First-principles characterization of ferromagnetic MnsGea for spintronic applications 
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In the active search for potentiaUy promising candidates for spintronic apphcations, we focus on 
the intermetallic ferromagnetic MnsGes compound and perform accurate first-principles FLAPW 
calculations within density functional theory. Through a careful investigation of the bulk electronic 
and magnetic structure, our results for the total magnetization, atomic magnetic moments, metallic 
conducting character and hyperfine fields are found to be in good agreement with experiments, 
and are elucidated in terms of a hybridization mechanism and exchange interaction. In order to 
assess the potential of this compound for spin-injection purposes, we calculate Fermi velocities and 
degree of spin-polarization; our results predict a rather high spin-injection efficiency in the diffusive 
regime along the hexagonal c-axis. Magneto-optical properties, such as L2,3 X-ray magnetic circular 
dichroism, are also reported and await comparison with experimental data. 
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INTRODUCTION 

Mn-doped Ge has recently been proposed as a promis- 
ing candidate in the challenging field of diluted magnetic 
semiconductors (DMS)0, which aims at combining 
information logic and storage. For example, epitaxial 
single crystal films of Mna;Gei_2: {x < 8-10%) grown on 
GaAs(OOl) and Ge were found to exhibit Curie tempera- 
tures over the range 25 to 116 K, combined with a p-type 
semiconducting behaviour Q. Many efforts are presently 
devoted toward increasing the transition temperature up 
to or above room temperature. Within this framework, 
one possible way is to increase the concentration of mag- 
netic impurities. However, one of the key issues in DMS 
is indeed the solubility of Mn in the semiconducting host: 
it is well known that beyond a certain critical Mn concen- 
tration (typically of the order of a few percent in III-V 
hosts), a tendency toward clustering and phase separa- 
tion occurs, thereby limiting the homogeneity and growth 
control that are strictly required for materials to be used 
in spintronic applications. This tendency was observed 
also during Mn-alloying of Ge samples: Mn^Gcj^ precipi- 
tates were detected during out-of-equilibrium growth]^. 

Intermetallic compounds of Mn and Ge occur in several 
different stoichiometrics and crystallographic phasesQ, 
most of which are antiferromagnetic or ferrimagnetic 
with rather low ordering temperatures. However, 
MusGca shows ferromagnetism with a Curie temperature 
of ~300 K, along with a uniaxial magnetic anisotropy 
along the c axis of the hexagonal crystal structure (see 
below) [el H Is! . Ferromagnetic MnsGea thin films grown 
epitaxially on Ge(lll) by means of solid-phase epitaxy ^S] 
exhibited metallic conductivity and strong ferromag- 
netism up to 296 K — thus holding out promise for use 



in spin injection. Moreover, very recently point contact 
Andreev reflection spectroscopy was used to measure the 
spin-polarization of MnsGea epilayers|9| and the results 
were compared with calculated values within the density 
functional theory. The discrepancy between the exper- 
imental and predicted spin-polarization was attributed 
to the extreme sensitivity of calculated results to the 
crystallographic structure, as well as to possible Mn de- 
ficiencies in Mn^Gcs samples. 9j Finally, it was shown 
experimentally, lOj upon C doping (with carbon intersti- 
tially incorporated into the voids of Mn octahedra of the 
MnsGea compound), that the Curie temperature, Tc, 
dramatically increased: MnsGcsCa; films for C concen- 
tration X > 0.5 showed Tc ^ 680 K. 

So far, very little is known theoretically about MnsGea 
; in particular, a careful investigation from first-principles 
of the magnetic interactions and chemical bonding be- 
tween Mn and metalloid atoms is still lacking. In this 
work, we perform a comprehensive study of MnsGca 
within density functional theory; in particular, in Sec^ 
we report the technicalities related to the structure and 
to the computational approach. The electronic structure, 
as well as the related magnetism, is discussed in Scc„ in 
terms of band structure, orbital and spin magnetic mo- 
ments, hyperfine fields and magnetic-circular dichroism 
spectra. Conclusions are drawn in SeCn 



STRUCTURAL AND COMPUTATIONAL 
DETAILS 

Our calculations were performed using one of the 
most accurate available density functional theory (DFT) 
methods, namely the all-electron full-potential linearized 
augmented plane wave (FLAPW) [ll| approach. The 



2 



generalized gradient approximation (GGA) according 
to the Perdew-Becke-Erzenhof scheme 12] was used for 
the exchange-correlation (XC) potential. This choice 
was suggested by the more accurate treatment of this 
exchange-correlation functional for magnetic compounds 
with respect to the local spin-density approximation 
(LSDA) [13] ; however, in order to test the resulting effects 
of a different (XC) parametrization and for the evalua- 
tion of the hyperfine fields, we also performed some calcu- 
lations using the von Barth-Hedin 14] functional within 
the local spin density approximation (LSDA). We used 
plane waves with wave vector up to K^ax = 3.8 a.u., lead- 
ing to about 1500 basis functions, whereas for the poten- 
tial and the charge density we used an angular momen- 
tum expansion with Imax < = 8. The Brillouin zone sam- 
pling was performed using 60 special fc-points in the irre- 
ducible wedge, according to the Monkhorst-Pack scheme 
(isf. The muffin tin radii, Rmt, for Mn and Ge were cho- 
sen equal to 2.37 a.u. and 2.0 a.u., respectively. In order 
to evaluate the effects of the orbital contribution to the 
magnetic moments, the calculations were performed with 
and without the spin-orbit coupling (SOC) included in 
the Hamiltonian|l6| . 

For the purpose of calculating the electronic group ve- 
locity v(k)= {l/fi)[de{'k)/dk], the eigenenergies e(k) over 
a set of 150 k points were used for a spline fitting of the 
bands over the Brillouin zone The resulting inter- 

polating Fourier series was then used to calculate the 
required energy derivative. A similar approach has been 
followed to calculate the electronic plasma frequency: 
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FIG. 1: (a) Perspective, (b) top and (c) side views of MnsGes. 
Black, white and grey spheres denote Mnl, Ge and Mn2 
atoms. The unit cell is also shown. 



where <> denotes the Fermi surface average. Given the 
hexagonal symmetry, the quantities, ujpxx — ^p\\ and o^pzz 
— tj^p±, will be evaluated. 

According to Forsyth and BrownQ, intermetallic 
MnsGea has an hexagonal crystal structure of DBs type 
(space group PGa/mcm), with experimental cell dimen- 
sions at room temperature a = 7.184 A and c = 5.053 
A. The atomic positions are: 



Mnl in 4(d) site: ±(i, 
Mn2 in 6(g) site: ±(x, 0, i; 0, a;, i 
0.2397 

Ge in 6(g) site: ±(a;, 0, \; 0, x, i; 

0.6030 
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Starting with the experimental equilibrium parame- 
ters, we checked that the calculated internal atomic forces 
were negligibly small and that the minimum total energy 
was obtained for the a and c value reported in Refl^. 
This confirmed that the FLAPW method as well as the 
GGA parametrization accurately reproduce the experi- 
mental structural properties for compounds with a high 
concentration of magnetic atoms. In Figure ^ we show 



the perspective, top and side views of the crystal. It is 
evident that there are two different atomic planes per- 
pendicular to the [001] direction: the first contains only 
Mnl atoms (at z = and z = c/2, equivalent by symme- 
try) forming an hexagonal two-dimensional lattice; the 
second contains Mn2 and Ge atoms (at 2; = c/4 and 
z = 3c/4, equivalent by symmetry). We recall that in 
the complex MnsGea structure, Mnl and Mn2 atoms 
have different coordinations; in particular0 , the nearest 
neighbours (nn) of each Mn atoms are arranged as: 

• Mnl has (i) two (six) Mnl (Mn2) nearest-neighbors 
at 2.522 (3.059) A and (m) six Ge at 2.534 A; 

• Mn2 has two Mn2, four Mn2 and four Mnl at 2.976, 
3.051 and 3.059 A, respectively and (ii) two Ge, 
one Ge and two Ge at 2.482, 2.606 and 2.762 A, 
respectively. 

In this configuration, the enthalpy of formation, A.Hf, 
of MnsGea is evaluated with respect to the stable phases 
of Mn ([001]-ordered antiferromagnetically fee) and Ge 
(in the zincblende phase). Our GGA calculated value. 
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FIG. 2: PDOS of (a) Ge, (b) Mnl and (c) Mn2 the total 
density of states is shown in panel (d). Majority (minority) 
spin components are shown in the positive (negative) y-axis. 
The Ef is set to zero of the energy scale. 



A Hf = 0.84 eV/formula-unit {i.e. ~0.1 eV/atom), 
shows that MnsGea is a quite stable compound. 



ELECTRONIC AND MAGNETIC PROPERTIES 



Density of states and band structure 



tivity experimentslg. The Ge atom shows slightly dif- 
ferent PDOS for up and down spins, consistent with the 
small negative magnetic moment (see below). In the en- 
ergy range considered, the largest contribution is due top 
states. On the Mn sites, as expected, contributions from 
s and p states (not shown) are negligible and the PDOS 
is essentially dominated by the 3d states. As a difference 
with the case of Mn impurities in Ge0], we here have 
a much smaller hybridization between Ge p and Mn d 
states. In fact, we find the Ge p states at higher binding 
energy while the Mn d dominate the region close to the 
Fermi level, showing a quite large Mn-Mn interaction. 

Differences in the PDOS of Mnl and Mn2 are partic- 
ularly evident in the majority spin occupied states; in 
particular, from the analysis of the peaks on different 
atoms located at the same energies, we can infer that: i) 
the Mnl features at ~-0.7 eV and in-between -2 and -3 
eV are due to Mnl-Mnl interactions; ii) the feature at 
~-l-2 eV is due to Mnl-Mn2 interactions; Hi) the Mn2 
feature at -3.5-2.5 eV can be ascribed to Mn2-Mn2 in- 
teractions. The high-binding energy range (<-3.5 eV) 
shows hybridization of both Mn atoms with Ge. Simi- 
larly, minority states for binding energies greater than 
1.3 eV show common features for Mnl, Mn2 and Ge 
atoms, whereas the feature at around -1-0.5 eV results 
from Mnl-Mn2 hybridization. The unoccupied states, of 
interest for the discussion of magneto-optical properties 
(see below), are largely due to the minority spin compo- 
nent and only show minor differences between Mnl and 
Mn2. 

The band structure for the majority and minority spins 
is shown in Figure El We remark that the levels around 
Ef (see Figure |21 (b) and (d)) are rather dispersed for 
the majority spin channel, whereas they are more local- 
ized in the minority spin component; moreover, a non- 
dispersed region is evident for higher binding energies in 
the majority spin band structure (in the energy range 
from -0.7 eV to -3.5 eV), as well as in the unoccupied 
minority spin band structure. The dispersion around Ep 
in the up-spin component confirms the strong hybridiza- 
tion between Mnl and Mn2 states, and between Mn d 
and p Ge states. Roughly speaking, in fact, as already 
pointed out for the DOS, the region at higher binding 
energies {i.e. in the energy range between -7 eV and -1.5 
eV) is basically due to a large contribution from Ge; on 
the other hand, the levels around Ep are basically due to 
both Mnl and Mn2 in the minority spin channel, whereas 
a p Ge contribution hybridized with Mn states is evident 
in the majority spin bands. 



The projected density of states (PDOS) of the three 
different atomic types forming the unit cell and the to- 
tal DOS are reported in Figure El As for its conduct- 
ing character, MnsGea shows strongly metallic behav- 
ior in both the minority and majority spin components: 
this is consistent with recently reported electrical resis- 



Magnetic moments, spin and charge density 

The calculated total magnetization and the magnetic 
moments of the different atomic species are compared 
with their corresponding experimental values and re- 
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FIG. 3: (a) Majority and (c) minority band structure along 
the main symmetry lines. Panels (b) and (d) show a blow-up 
around the Fermi level. The Ef is set to zero of the energy 
scale. 



TABLE I: Magnetic moments within the atomic spheres and 
total magnetization per unit cell (in Bohr magnetons). The 
first two hues show the spin magnetic moments {ij,nosoc) ^ 
obtained from calculations without SOC within LSDA and 
GGA; the third, fourth and fifth lines show the spin {^J.soc)^ 
orbital (/^soc) ^'id total (fJ-'soc) magnetic moments as ob- 
tained including the SOC self-consistently within GGA. The 
magnetic moments of the atomic species are calculated within 
their muffin-tin spheres, whereas the total contribution in- 
cludes the contribution from the interstitial region. Exper- 
imental values (fiixp) from RefsQ, are shown in the last 
row. 





Ge 


Mnl 


Mn2 


Total 


^^''Nosoc (LSDA) 


-0.09 


2.14 


3.11 


26.7 


f^NOSoc (GGA) 


-0.11 


2.28 


3.22 


27.5 


Msoc (GGA) 


-0.11 


2.07 


3.12 




fi'soc (GGA) 





0.05 


0.035 




f/soc (GGA) 


-0.11 


2.12 


3.16 


25.9 


t 

f-^exp 




1.96" 


3.23" 


26'' 



a. Refl| 

b. Refif 



ported in TableHJ We show the magnetic moments (both 
spin and orbital contributions) with and without the in- 
clusion of SOC. As experimentally reported^, all the 
moments are ferromagnetically oriented along the crys- 
tallographic c axis. A comparison between LSDA and 
GGA does not show any significant differences as far 
as the general magnetization distribution over different 
atomic sites is concerned; however, LSDA predicts lower 
atomic magnetic moments (by about 0.1-0.15 fis) than 
GGA, resulting in a LSDA total magnetic moment per 
unit cell that differs by almost one Bohr magneton with 
respect to its GGA counterpart. 

According to previous neutron scattering 
experiments [12|J, the magnetic structure of MnsGea 
is found to reveal two Mn sublattices (Mnl in a four- 
fold and Mn2 in a six-fold position) with different 
magnetic moments. It was suggested that Mn2 carries 
the larger moment, in agreement with zero-field NMR 
measurements |20j . Our results are in excellent quanti- 
tative and qualitative agreement with these results. As 
suggested in Ref|3, the lower Mnl magnetic moment 
is due to the different Mn coordination and to direct 
Mn-Mn interactions at a rather short distance (recall, 
in fact, that every Mnl atom in MnsGes has 2 Mnl at 
a distance of 2.52 A): for a large set of Mn intermetallic 
compounds, an analysis of the Mn magnetic moment 
vs the nearest neighbor distance^ showed that below 
a "threshold" separation of 3.1 A, a moment reduction 
of ~ 2 /is/ A per Mn neighbor occurred with respect 
to the atomic value (5 /^s) of the Mn^+ ion, therefore 
accounting for the small magnetic moment at the Mnl 
site. 

It is remarkable that the inclusion of the orbital mo- 
ments (which are not completely negligible on both the 
Mnl and Mn2 sites) largely improves the agreement with 
experiment: the total magnetic moment is in excellent 
agreement with the experimental saturation magnetiza- 
tion as obtained by Kappelj2lt . Moreover, it was sug- 
gested by Forsyth and Brownj^ that some spatially dif- 
fuse reverse magnetization (in contrast to localized mo- 
ments centered on the Mn atoms) exists in MnsGea. This 
is consistent with our calculated negative moment on Ge 
sites. As expected, the Mn magnetic moments are due to 
d states, whereas the small induced moment on Ge is due 
to p state polarization. In order to further investigate this 
issue and to better show the bonding, we plot the charge 
and spin density on two different planes perpendicular to 
the [0001] axis. The charge density in Figure0](a) shows 
fairly "isolated" Mnl atoms; a somewhat stronger inter- 
action occurs between Ge and Mn2 atoms, as shown by 
the presence of charge in the bonding regions in Figure 

EKc). 

As regards deviations from spherical symmetry in the 
spin density, we remark that Mn2 shows an almost spher- 
ical shape, whereas the spin density around Mnl is ex- 
tended towards the six surrounding Ge ligands. It is ev- 
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FIG. 4: (a) Valence charge density and (b) spin density in 
a plane perpendicular to the [0001] direction and containing 
Mnl atoms, (c) Valence charge density and (d) spin density 
in a plane perpendicular to the [0001] direction and contain- 
ing Mn2 and Ge atoms. Shading for the atomic spheres are 
consistent with Figure 1: white, black and grey circles denote 
Ge, Mnl and Mn2 atoms, respectively. In panels (b) and (d) 
solid (dashed) lines show the positive (negative) contribution 
to the spin density. 



idcnt from Figure ^d) that a negative spin-density sur- 
rounds the metalloid Ge atoms, consistent with the neg- 
ative magnetic moment reported in Table ^ Moreover, 
the diffuse negative spin density shown in Figure 01 (b) 
and (d) confirms the experimental results. 

Hyperfine fields 

We focus next on hyperfine fields and compare our 
predicted values with experimental dataj^. As is well- 
known, the hyperfine field[23 of an atom is the magnetic 
field at the atomic nuclear site produced by the electrons 
in the solid and can be probed using Mossbauer spec- 
troscopy or nuclear magnetic resonance to provide valu- 
able information on the electronic and magnetic prop- 
erties of the compound. It consists of several contribu- 
tions: (i) the leading term due to the Fermi contact in- 
teraction which is proportional to the spin density at the 
nucleus [2^. 

Kf^ln^iipm-pim, (2) 

in the scalar relativistic limit; (ii) an orbital term which, 
according to Abragam and Pryce|2^ can be expressed 
as: 

iir/ - 2^.3 < >i (3) 

where < r^^ >i is the average expectation value of 
of the radial wave function and /i' is the orbital mag- 
netic moment; (iii) a dipolar term Hf^^ . Whereas only 
s electrons contribute to the Fermi contact term, elec- 
tronic states with I 7^0 contribute to the latter terms. In 
scalar or non-relativistic calculations, the orbital angu- 
lar momentum is quenched and Hf^^jf'^O; however, when 
spin-orbit is included, this term can be non-vanishing. 
In our case, the p {I = 1) contribution to the orbital mo- 
ment is negligible (<10^'^ /is) and the orbital term will 
be therefore evaluated only for I — 2 {d states) . The dipo- 
lar contribution is normally small in bulk systems and is 
therefore neglected. 

In Table im we report our calculated values for the core 
and valence contributions to the Fermi contact hyperfine 



TABLE II: GGA calculated Fermi contact hyperfine fields 
broken down into the core (JY^'^"'^) and valence (-ff^^'""') 
contributions (values in parenthesis denote LSDA calculated 
values). The total Fermi contact hyperfine fields (-ff)J^'*°'), the 
average < > for the d states and the orbital hyperfine 
field (-ffh/*) ^re also shown. The total hyperfine field Hj^j = 
H'f^j^°^ + -f^h/'' is compared with available experimental data 
(magnitude of Jj"f^^')f20|. Values of hyperfine fields (< >) 
are expressed in kOe {a^^ , where ao is the Bohr radius). 





-rrct .core. 


Trct.val 


rjCt.tot 


< r ^ > 


TTorb 






Mnl 


-316 
(-279) 


111 

(76) 


-205 
(-203) 


2.14 


12.9 


-192 


195 


Mn2 


-459 
(-417) 


107 
(70) 


-352 
(-347) 


2.18 


9.5 


-342 


399 



field, along with the orbital and total contributions, com- 
pared with experimental values. The dominant exchange 
polarization of the core electrons has to be taken into 
account, showing the need for an all-electron method[Tl| 
when dealing with hyperfine fields. As discussed in previ- 
ous theoretical work for transition metals. [2^ in Mn5Ge3 
the separation of the negative core and positive valence 
contributions to the Fermi contact hyperfine field high- 
lights these two opposite terms (cf. Table Hill . The large 
negative core contribution can be attributed to the at- 
traction of the majority spin electrons towards the spatial 
region of the spin-polarized d shell |22j which produces 
the excess of minority spin electrons at the nucleus. In 
order to evaluate the effects of a different parametriza- 
tion for the exchange-correlation potential, we compare 
the Fermi contact term within LSDA and GGA. The to- 
tal (core-|-valence) contribution is very similar; however, 
both the separate core and valence terms have a larger 
magnitude within GGA. Note that the core polarization 
per unit spin moment for Mn 1 and Mn 2 in both LDA 
and GGA is constant - as is expected from the exchange 
polarization mechanism but with a somewhat different 
constant, namely ~ 130 kG//iB and ~ 140 kG//iB within 
LSDA and GGA, respectively. These values are pretty 
similar to the values obtained for the Mn-based Heusler 
compounds 24] (-- 140 kG/fis within LSDA and ~ 150 
kG/nB within GGA). 

As far as the orbital contribution is concerned, we 
point out that this term is, as expected, much smaller 
than the Fermi-contact term and very small, because 
the unquenched orbital moment is so small, as is usual 
for Mn. Its magnitude could be slightly underestimated 
due to the well-known failure of SDFT in determining or- 
bital magnetic moment; however, this error is expected 
not to dramatically change the final value of the total hy- 
perfine field. In particular, we point out that the inclu- 
sion of the (positive) orbital term improves (worsens) the 
agreement with experiment in the case of Mnl (Mn2). In 
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fact, by means of zero-field NMR and specific heat mea- 
surements, the magnitude of the experimental effective 
nuclear fields were determined as f95 kOe at the 4(d) 
Mn site and 399 kOe at the 6(g) Mn sitefl^. The agree- 
ment of the calculated values {Hlf{Mni) = -192 kOe 
and 7J^y*(Mn2) = -342 kOe) with experiments is seen to 
be reasonably good. 

FERMI VELOCITIES AND DEGREE OF 
SPIN-POLARIZATION 

Since MnsGca has been suggested as a potential spin- 
injector, it is useful for device applications to give in- 
formation about transport properties, in terms of Fermi 
velocities and spin-polarization. 

Recall that various definitions of spin-polarization P 
have be en p roposed, each of them to be used in different 
regimes [2^. The most natural and popular definition 
involves the DOS at Ep and is probed, for example, in 
spin-polarized photoemission measurements: 

However, in transport measurements (see for example 
Andreev reflection 26j ) , Fermi velocities are of course rel- 
evant quantities and should therefore be involved in the 
spin-polarization definition 25] . In particular, for low re- 
sistance ballistic contacts, the appropriate definition is: 

Pi = (< N{Ef)-vf >t - 

< N{Ef) ■ VF >i)/(< N{Ef) • >t + 

< N{Ef) ■ VF >i) 

whereas for large barrier and/or diffusive current the cor- 
rect definition is: 

P2 = {<N{Ef)-vI >t - 

< N{Ef) ■ vl >i)l{< N{Ef) ■vl>^ + 

<N{Ef)-vI >^) 

where <> denotes the Fermi surface average. 

In Figure|31we show our GGA-calculated in-plane (per- 
pendicular to the c axis) and out-of-plane (parallel to the 
c axis) velocities as a function of energy. 

Moreover, in the upper part of Table IIIII the corre- 
sponding quantities evaluated at the are reported, 
along with plasma frequencies. According to these val- 
ues, we have calculated the Pqj Pi and P2 values (the lat- 
ter two for different in-plane and out-of-plane directions) 
reported in Table ITTll (lower part). Due to the numeri- 
cal uncertainties (related to the k-point sampling, wave 
function cut-offs, etc.), we estimate an error on the spin- 
polarization of ±5-10 %. Within this error, our values 
are consistent with similar values recently obtained us- 
ing a different DFT method and slightly different lattice 
parameters, js] 

As pointed out in the discussion of the electronic prop- 
erties, the Mn heavy d bands in the majority channel are 
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FIG. 5: In-plane (solid) and out-of-plane (dashed) velocities 
as a function of energy for majority and minority spin elec- 
trons (positive and negative y axis, respectively). 



TABLE III: Upper part: GGA-calculated relevant quantities 
at E_F for up and down spin-channels (first and second line, re- 
spectively): Density of states {N{Ef), in states/eV), in-plane 
velocity {vf\\, in 10® cm/sec), out-of-plane velocity («f±, in 
10® cm/sec), in-plane plasma frequency {iOp^\, in eV) and out- 
of-plane plasma frequency (tJpx, in eV). Spin-orbit coupling is 
not included. Lower part: Degree of spin-polarization (DSP) 
as measured in photoemission measurements (-Po), ballistic 
transport (in-plane and out-of-plane, Piy and Pi±, respec- 
tively) and diffusive transport (in-plane and out-of-plane (P211 
and P2±, respectively.) - see text for definitions. For compar- 
ison, we also show similar quantities obtained for ferromag- 
netic hep Co. 





Compound 


N{Ef) 


vf\\ 


Vf± 




UJp± 


Spin-up 


MusCea 


3.3 


19.3 


22.6 


2.0 


2.4 


Spin-down 




7.9 


6.7 


6.4 


1.1 


1.1 


Spin-up 


Co 


0.3 


47.9 


38.2 


5.0 


4.0 


Spin-down 




1.5 


14.0 


16.4 


3.3 


3.8 




Compound 


Po 


^111 


Pix 


P211 


P2± 


DSP 


MnsGea 


-41% 


8% 


18% 


54% 


67% 




Co 


-67% 


-19% 


-37% 


40% 


4% 



almost fully occupied and the Fermi level also crosses Ge 
p states, which are light states with an appreciable veloc- 
ity. On the other hand, in the minority spin-channels, the 
DOS shows contributions from both Mn heavy d and Ge 
p states. As clearly shown in Table Hill the Fermi DOS 
in the minority spin-channel is larger - by a factor of ^^2 
- than in the majority spin-channel, whereas the Fermi 
velocities (both in-plane and out-of-plane) are larger by 
a factor of ^^3 for majority spins. This leads to a nega- 
tive spin-polarization Pqj but to a positive current (see 
positive values of Pi), as also noted in 3d metals. [j^j As 
for the plasma frequencies, we point out that there are 
only slight differences between in-plane and out-of-plane 



calculated values, whereas larger differences between mi- 
nority and minority spins emerge: uij, is about two times 
larger than w^. 

As shown in Figure [S] the anisotropy of the in-plane 
with respect to the out-of~plane velocity is quite evi- 
dent: although the behavior as a function of energy is 
overall similar, they slightly differ in the case of major- 
ity spins, at about -leV and in the relevant energy range 
around the Fermi level (from -0.2 to 0.5 eV in the ma- 
jority channel). Moreover, as shown in Tabic ITTll the 
different definitions of spin polarizations result in largely 
differing values. Both the anisotropy as well as the dif- 
ferences among Pq, Pi and P2 should help in exploiting 
this compound for spin-injections purposes in the most 
appropriate transport regime and along the most favor- 
able growth direction. In particular, we remark that with 
P2± ^ 70%, most of the current along the c axis in the 
diffusive (Ohmic) regime is therefore carried by major- 
ity spins. Of course, this picture might be modified in 
the presence of a junction (such as MnsGea/Ge), where 
interface states can modify the electronic structure and 
velocities with respect to the ideal bulk situation consid- 
ered here. 

Finally, it is useful to compare calculated spin-injection 
efficiencies and Fermi velocities for MnsGea with those 
obtained for ferromagnetic hep Co[2_7|, a widely used ma- 
terial in high density magnetic recording media. Our re- 
sults are shown in Table IIIII Interestingly, the cobalt 
Fermi velocities (and related plasma frequencies) are 
larger by a factor of 2-3 compared to MnsGca , also 
showing a larger anisotropy. As far as the Co DSP is 
concerned, we remark that, similarly to MnsGca , Pq is 
negative; however, Pi - both in-plane and out-of-plane 
- is also negative, at variance with MnsGea. As a fi- 
nal remark, we observe that the Co spin polarization in 
the ballistic regime is slightly higher than in MnsGca , 
whereas it is definitely lower in the diffusive regime, hold- 
ing promise for MnsGea as an efficient spin-injector. 



X-RAY ABSORPTION AND MAGNETIC 
CIRCULAR DICHROISM 

The formalism within band theory to calculate the 
cross section for the absorption of incident light is dis- 
cussed in Ref l2^ and is briefly reviewed here. For dipole- 
excited transitions, the cross section can be obtained as: 

an{E) = [ I < <fc\p\^, > {""SiE - {E, - E,))dk (4) 

where n = z, ± represents the photon polarization (i.e. 
incident light polarized vertically along the direction of 
magnetization (z) or left- (-I-) or right-circularly (— ) po- 
larized), p is the momentum operator, and {Ec and 
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FIG. 6: (a) L2,3 X-ray absorption; (b) L3 X-ray MCD; and (c) 
L2 X-ray MCD. The insets in panels (b) and (c) show MCD 
spectra convoluted with an energy broadening of 0.7 eV. In 
all panels, the solid (dashed) line shows the contribution of 
the Mnl (Mn2) site. 



Ey) denote initial core and final valence states (eigenener- 
gies), respectively. Therefore, the X-ray magnetic circu- 
lar dichroism (XMCD) can be obtained as dm = (7+ — (t_, 
whereas the X-ray absorption spectrum (XAS) is calcu- 
lated as at = <y+ ~\- a - -V ((Tq denotes the absorption 
cross section for incident light polarized vertically along 
the direction of magnetization, here chosen as coincident 
with the z hexagonal axis). Starting from the converged 
FLAPW ground state, the spin-orbit coupling (SOC) was 
treated in a second variational way .16] to obtain the XAS 
and XMCD spectra, using up to 432 k-points in the full 
Brillouin zone. We used a 0.25 eV Lorentzian broadening 
to smooth the calculated spectra, in order to take lifetime 
effects into account (see below). 

The energy dependence of L2.3 fm and at are shown 
in Figure O for the two Mn atomic types. As is usual for 
3(i metals, j28j and pointed out above in the PDOS dis- 
cussion, most of the majority spin bands are located be- 
low Ep, so that photon induced transitions occur mainly 
to the unoccupied minority spin-bands. We haven't in- 
cluded any self-energy correction, so we expect the cal- 
culated binding energy of the 2p states to be underesti- 
mated by several tenths of an eV. The energy difference 
between the L2 and is edge represents the size of the 
spin-orbit splitting of the 2p core states, and is usually 
found to be in good agreement with experiment. For 
MnsGea it is estimated to be ^^Joc ^ ^0.4 eV. This 
value is similar to other theoretical results for Mn based 
alloys, such as the Heusler PtMnSb and NiMnSbjl^, as 
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well as to other experimental data obtained for DMS|30(. 
As expected from the quite low symmetry of the hexago- 
nal lattice (and the related difference of the out~-of~plane 
z direction compared to the in-plane x, y direction), ctq is 
quite different from ^(cr+ + ct_) (not given in Figure O. 
The XAS spectrum shows tails extending to high energy 
(>10-15 eV) with respect to the main absorption peaks, 
whereas the XMCD becomes almost negligible at 5-7 eV 
above the absorption edge. For the energy broadening 
value used (0.25 eV), both the absorption and dichro- 
ism spectra show a quite rich structure; in particular, 
the XAS and - even more markedly - the XMCD spec- 
tra show a double peak structure (related to the peculiar 
features in the unoccupied density of states, see Fig|21 
(b) and (c)), followed by a smaller bump at about 3.5 
eV above the absorption edge. As far as the compari- 
son between Mnl and Mn2 is concerned, we remark that 
the XAS spectrum shows similar features, whereas the 
XMCD shows different amplitudes of the second peak, 
which is more marked in Mn2 compared to Mnl. This 
difference in the XMCD intensity is consistent with the 
larger Mn2 magnetic moment compared to Mnl. Since 
the broadening used (0.25 eV) is lower than the common 
experimental resolution, we also show in the insets of 
FigEl (b) and (c) the same XMCD spectra obtained us- 
ing a Lorentzian broadening equal to 0.7 eV; in this case, 
the double-peak feature cannot be resolved, whereas the 
larger amplitudes observed for Mn2 compared to Mnl in 
the high energy range is still evident. 

To further investigate this issue, we recall that some 
important magneto-optical sum rules have been derived 
in recent years, which relate the integrated signals over 
the spin-orbit split core edges of the unpolarized XAS 
and of circular dichroism to ground-state orbital and spin 
magnetic moments |3]1 [33- .34] . The orbital and spin 
sum rules are expressed as: 



2/™,7Vft, 



< s. >= — 7 < r, > 



(5) 
(6) 
(7) 
(8) 
(9) 



where Nh is the number of holes in the d band, = 10 — 
n3d (with n^d determined by the d projected density of 
states inside each atomic sphere) . is the z component 
of the magnetic dipole operator: 

T, = 1/2 [a - 3 f (f • (t)1 



TABLE IV: Orbital and spin magnetic moments as deter- 
mined from sum-rules (SR) and self-consistently (SC) for 
Mnl and Mn2 atoms along with their ratio. 













SR 


SC 


SR 


SC 


SR 


SC 


Mnl 


0.03 


0.05 


1.9 


2.07 


0.015 


0.024 


Mn2 


0.02 


0.035 


2.4 


3.16 


0.01 


0.01 



with (T denoting the vector of Pauli matrices and is re- 
lated to the non-spherical charge and spin density. is 
vanishing for cubic systems, whereas it is not necessarily 
negligible in hexagonal systems. However, preliminary 
calculations suggest that is pretty small (at most of 
the order of 0.03) and is therefore neglected here. The 
orbital and spin magnetic moments are /i' = —fiB < Iz > 
and /i* = —^J■B < Sz >, respectively. Their ratio, as de- 
rived from sum rules, reads as: 



<lz> 
<Sz> 



3 Is 
2I„, 



7<T,>It 



1 -1 



2 /,„ Nh 



(10) 



Experimentally, since < > is very hard to access, 
it is generally neglected when considering sum rules; in 
this case the experimental uncertainties - related to fixing 
somehow the number of holes or to calculating /( - 
drop out. It has therefore been suggested 0| that an 
accurate estimate of the /i'/M* ratio can be obtained from 
the 2/m/3/s ratio. 

In Table IIVI we report our calculated values for the 
orbital and spin magnetic moments, as well as their ra- 
tio, as determined from sum rules (SR) for Mnl and 
Mn2 atoms. For comparison, we also report the same 
values as self-consistently calculated (cfr. Tabled). As 
expected, the general trends of spin and magnetic mo- 
ments for Mnl and Mn2 atoms are qualitatively well 
reproduced. However, quantitatively, the SR magnetic 
moments are generally underestimated with respect to 
the self-consistently calculated values (ranging from ~10 
percent in the case of spin magnetic moments up to 40 % 
in the case of the small orbital moments). This has been 
ascribed |33| to the several approximations made in deriv- 
ing the sum rules [3^, among which the most serious are 
i) ignoring the interatomic hybridization, ii) neglecting 
the p ^ s transitions 3^, and in) ignoring the exchange 
splitting of core levels. Therefore, first-principles calcula- 
tions of both XMCD spectra and ground-state magnetic 
moments are crucial for a careful study of the MnsGca 
compound and, eventually, for a quantitative interpreta- 
tion of future experimental results. 



SUMMARY 

In the search for new compounds to be used for efficient 
spin-injection in spintronic devices and following the re- 
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cent suggestion of Mn5Ge3/Ge(lll) system as a promis- 
ing system, we have presented a careful first-principles 
FLAPW investigation of the electronic, transport, mag- 
netic and magneto-optical properties of bulk MusGes. 
Our results show that the two Mn sites have different 
magnetic moments, leading to a total magnetization of 
26 /is per unit cell; the conducting character is strongly 
metallic with states around the Fermi level essentially 
due to Mn-Mn interactions. Our theoretical predictions 
were carefully analyzed in terms of the underlying elec- 
tronic and magnetic structure and shown to be in excel- 
lent quantitative and qualitative agreement with avail- 
able experimental results. The most favorable condition 
for spin-injection purposes is predicted to be in the dif- 
fusive regime along the hexagonal c axis, where a rather 
high spin-polarization is obtained. 

We thank Prof. Ruqian Wu and Prof. Sandro Mas- 
sidda for assistance provided with the calculation of the 
magnetic dipole operator and spline fitting procedure, 
respectively. Useful discussions with Dr. Steven C. Er- 
win are gratefully acknowledged. Work in L'Aquila sup- 
ported by INFM through Iniziativa Trasversale Calcolo 
Parallel© and PAIS-GEMASE project. 
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